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Abstract 

We investigate the non-minimal couplings between the electromagnetic fields and gravity through 
the natural logarithm of the curvature scalar. After we give the Lagrangian formulation of the 
non-minimally coupled theory, we derive field equations by a first order variational principle using 
the method of Lagrange multipliers. We look at static, spherically symmetric solutions that are 
asymptotically flat. We discuss the nature of horizons for some candidate black hole solutions 
according to various values of the parameters Rq and a%. 
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I. INTRODUCTION 



We consider the non-minimal couplings of gravitational and electromagnetic fields de- 
scribed by a Lagrangian density. Such couplings may occur near compact astrophysical 
objects which has high mass density such as the neutron stars or black holes. The non- 
minimally coupled electromagnetic fields to gravity in RF 2 form were extended and classified 
[H [2] to gain more insight into the relationship between space-time curvature and charge con- 
servation. They were also obtained from a calculation in QED of the photon effective action 
from 1-loop vacuum polarization on a curved background [3] and from Kaluza-Klein reduc- 
tion of a five- dimensional i? 2 -Lagrangian [U |5]. A non-minimally coupled Einstein-Maxwell 
Lagrangian in general may involve field equations of order higher than two. However, the 
type of non-minimal couplings we consider are obtained by the reduction of the Euler- 
Poincare Lagrangian in five dimensions to four dimensions and involve at most second order 
field equations [6J. A three parameter family of non-minimally coupled Einstein-Maxwell 
field equations was studied in [7]. 

Recently, in the context of primordial magnetic fields present during the reheating epoch 
of the universe, the i? m F 2 -type couplings were discussed in [8]. The modified gravity with 
InR terms [9] has also got much attention, since it could explain the observed acceleration 
of the universe and the dark energy concept without using any exotic fields. The behavior 
of the rotational velocities of test particles gravitating around galaxies is investigated in 
modified gravity in the context of non-minimal matter couplings [TU] ITT] . In this context, 
the modified /(i?)-Maxwell gravity with I(R)F 2 coupling terms [121 [13] an d /(G)-Maxwell 
gravity in which non-minimal coupling between electromagnetic fields and a function of 
Gauss-Bonnet invariant [H] were proposed to explain late-time cosmic acceleration. There 
are many other theories of gravity with non-minimal couplings [15], [16]. However, the 
non-minimal couplings with electromagnetic fields are not investigated in sufficient detail. 
Especially, finding spherically symmetric solutions is not an easy task for such theories 
[IT], [18]. Furthermore, any arbitrary non-minimal coupling may not give rise to solutions 
satisfying physical asymptotic conditions and observations in solar and cosmological scales. 
Therefore, here we propose a non-minimal theory with a special InR coupling term. 

We first discuss a non-minimal Y(R)F 2 coupled Einstein-Maxwell theory using the al- 
gebra of exterior differential forms. We derive the field equations by a first order varia- 
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tional principle using the method of Lagrange multipliers. We choose in particular Y(R) = 
l-aiinR/Ra an d ^ or static, spherically symmetric solutions that are asymptotically flat. 
In our model, the Y(R)F 2 coupled term in the Lagrangian leads to modifications both in the 
Maxwell and Einstein field equations. The modifications in the Maxwell equations can be 
related with the polarization and the magnetization in a specific medium. The non-minimal 
couplings also give important modifications to the structure of a charged black hole. These 
may shed light on some important problems of gravity such as dark matter and dark energy 
without introducing a cosmological constant or any other type of scalar fields. 



II. FIELD EQUATIONS OF THE NON-MINIMALLY COUPLED THEORY 



We will derive our field equations by a variational principle from an action 

I[e a ,oo\F]= [ L = [ (1) 

where {e a } and {oo a b } are the fundamental gravitational field variables and F is the electro- 
magnetic field 2-form. The space-time metric g = r] ab e a <8> e b with signature ( — h +- 1-) and 
we fix the orientation by setting *1 = e° A e 1 A e 2 A e 3 . Torsion 2-forms T a and curvature 
2-forms R a b of spacetime are found from the Cartan-Maurer structure equations 

T a = de a + u\ A e\ (2) 

R\ = dw\ + u a c A u\. (3) 
We consider the following Lagrangian density 4-form; 

L=±R*1-±Y(R)FA*F, (4) 

where n 2 = 8nG is Newton's universal gravitational constant (c = 1) and R is the curvature 
scalar which can be found by applying interior product i a twice to the curvature tensor R ab 
2-form. We use the shorthand notation e a A e b A • • • = e ab "' , and t a F = F a , i ba F = F ab , 
i a R a b = Rb, tb a R ab = R- The field equations are obtained by considering the independent 
variations of the action with respect to {e a }, {uJ a b } and {F}. The electromagnetic field 
components are read from the expansion F = \F ab e a A e b . We will confine ourselves to 
the unique metric-compatible, torsion-free Levi-Civita connection. We impose this choice 
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of connection through the constrained variations of the action by the method of Lagrange 
multipliers. That is, we add to the above Lagrangian density the following constraint terms: 

L c = (de a + cu a b A e b ) A A a + dF A //, (5) 

where A a 's are Lagrange multiplier 2-forms whose variation imposes the zero-torsion con- 
straint T a = 0. We also use a first order variational principle for the electromagnetic field 
2-form F for which the homogeneous field equation dF = is imposed by the variation of 
the Lagrange multiplier 2-form /i. 

The infinitesimal variations of the total Lagrangian density L + Lc (modulo a closed form) 
are given by 

L + L c = ^- 2 e a A R bc A *e abc + e a A ^Y(R)(t a F A *F - F A i a * F) + e a A DX a 

+e a A Y R (t a R b )(t b F A *F + F A t b * F) + ^u ab A (e b A \ a — e a A A 6 ) 

u ab A S ab - F A Y(R) * F + A a A T a - F A dfi. (6) 

where Yr — and the angular momentum tensor 

S a6 = --D[Y R (F ab * F + F b A i a * F - F a A i b * F - F A i ab * F)}. (7) 

The Lagrange multiplier 2-forms A a are solved uniquely from the connection variation equa- 
tions 

e a A \ b - e b A \ a = 2S afe , (8) 

by applying the interior product operator twice as 

A a = 2i b T, ba + -i bc Y, cb A e a . (9) 

We substitute the A a 's into the e a equations and after some simplifications we find the 
Einstein field equations for the extended theory as 

^R bc A *e abc + \Y{i a F A *F - F A i a * F)+Y R (i a R b )i b (F A *F) 

+ ^D[L b D(Y R F mn F mn )]A*e ab = 0, (10) 

while the Maxwell's equations read 

dF = , d*(YF) = 0. (11) 
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In terms of a local inertial coordinate system (x^), the master equations (10) and (11) turn 
out to be equivalent to the following equations that were obtained by Bamba and Odintsov 
[13J by a variation procedure with respect to the metric: 

where the electromagnetic energy-momentum tensor components 

T, u =Y(g^F^F ua - l -g^ v F a ^) 

+^{Y R F a ^R f , v + g^'VyV^RFapF ^ - V „V ' v [Y R F a pF aP }} , (13) 

and the Maxwell equations 

' --d^V^gPn = , ^^[v^FFH = 0. (14) 



V~9 V~9 
The effects of non-minimal couplings of the electromagnetic fields to gravity can be expressed 
by a constitutive tensor. For this purpose, Maxwell's equations satisfied by electromagnetic 
field F in an arbitrary medium can be written as 

dF = , *d * G = (15) 

where G is called the excitation 2-form. We have taken the source-free Maxwell equations. 
The effects of gravitation and electromagnetism on matter are described by G. We can 
complete this system writing the following linear constitutive relation 

G = Z(F) (16) 

where Z is a type-(2,2)-constitutive tensor. For the above theory, we have 

G = Y(R)F. (17) 

We can further introduce the polarization 1-form p = d— e — (Y — 1)ljjF and magnetization 
1-form m=b — h=(l — Y)lu * F relative to a time-like unit velocity vector field U of an 
inertial observer ( For more details see [T9l - f2~l~] ). 
Now we consider the following function Y(R) as a special case: 



1 

1 — a% In 



Y (R) = * I 5 ( 18 ) 



Ro 
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where a\ is a dimensionless coupling constant and Rq is a constant with the same dimension 
as that of R. We take a\ — to get back to the minimal Einstein-Maxwell theory. Further- 
more, if we look at the limit R — > Rq we again obtain the minimal Einstein-Maxwell case. 
It is interesting to note that in the limits as R goes either to zero or to infinity, the function 
Y(R) — y 0. That is, the effects of gravity become important and the electromagnetic effects 
can be neglected in these limits. On the other hand, as R approaches Rqc 1 ^ 1 , Y(R) becomes 
very large so that the electromagnetic effects dominate. Finally, the function Y(R) can be 
expanded as a power series in InR for a x 7^ and Icq In < 1 as 

1 00 p 

, flh H -E(^|)" (19) 

In this form the model resembles the RG improved theory proposed by Bamba and Odintsov 
in Ref.[13]. Following this article, we try to relate the asymptotic freedom in a non-Abelian 
SU(2) gauge theory with a non-minimal Maxwell-modified gravity by setting 



^=E(«ii*ir) n ( 2 °) 



where t is a renormalization-group parameter, g(t) is the running SU(2) gauge coupling 
constant and g = g(0). Therefore, if |a,i In 1 << 1, we obtain a similar kind of the RG 
parameter which has been proposed in 



III. STATIC, SPHERICALLY SYMMETRIC SOLUTIONS 

We seek static, spherically symmetric solutions to the field equations which are given by the 
metric 

g = -firfdt 2 + f(r)- 2 dr 2 + r 2 d6 2 + r 2 sm(6) 2 d(j) 2 (21) 
and a static electric potential 1-form A = V(r)dt. Then 

F = dA = V'dr A dt = Edr A dt. (22) 



A. Reduced Field Equations 



After a lengthy calculation we reduce the non-minimally coupled Einstein-Maxwell equa- 
tions (10) and (11) for the metric (21) and the electromagnetic 2-form (22) to the following 



system of equations: 

1 f 2 ' f 2 - 1 f 2 " f 2 ' 1 2 

- 2 (— + I — 2 -) - YrE 2 ^ + — ) + o YE " ~ [(& Y *Yf\'f ~ -f( E2Y R)' = o, 

k z r r z 2 r 2 r 
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1 f 2 ' f 2 - 1 f 2 f 2 ' 1 f 2 ' 2f 

- 2 {— + J — 2 —) ~ Y R E 2 { J — + J —) + -YE 2 - (E 2 Y R )'( J — + J-) = 0, (23) 
k z r r z 2 r 2 2 r 

i J?2'' /-2' i-2' f 2 _ 1 1 f 2 ' f 2 

^( V + ~ ] - YrE2{ ~ + — ] - 2 YE2 ~ K^'M - (^)'(t + y) = °> 

= (24) 

Here the curvature scalar 

^-/ 2// -V-^(/ 2 -l)- (25) 
q is the electric charge determined by the Gauss integral 

-*-/" *G = —[ Y(R)E(r)r 2 sin Odd A d<p = q. (26) 

B. Exact Solution 

At this point, we will make a simplifying assumption and consider solutions that satisfy 

R = S (27) 



where C is a constant to be fixed. Then, the field equations (23) and (24) simplify to 



R, a-[K 2 E 2 , . , 

2 ' ^l-.,.. (29) 

' _ 2 



£r 2 = g. (30) 



Thus, for ai 7^ 0, we can solve the unknown functions for the non-minimally extended 

Einstein-Maxwell theory with C = a\n 2 q 2 and find 

2M aiK 2 q 2 , r «V(l + 5ai) 
— + In - + \ 9 



f(r) = 1 - — + ^ In - + (31) 



4oai In — 

*W = £ + -Hr- 31 - (32) 

ro is an integration constant satisfying the relation r$ = ttl ^ g . The above solution is 
asymptotically flat. It is interesting to note that the electric field changes sign (due to 
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polarization effects) at r = roe 4a i that corresponds to the value R = Roe 1 ^ 1 . On the 
other hand, the displacement vector field is an inverse square force field and hence does 
not change sign. We observe that both the curvature scalar and the quadratic curvature 
invariant *(R a b A *R ab ) for our solutiuon are singular at the origin r = 0. 



C. Horizons and Asymptotic Behavior 

The roots of the metric function f 2 (r) are determined by the intersection points of a 
logarithmic curve and a parabola. In order to find these points we look at the numerator of 
the metric function: 

2 2 

n{r) = r 2 - 2Mr + —(1 + 5ai) + a X K 2 q 2 ln— = r 2 f 2 {r). (33) 

4 ro 

The derivative of the function or n'(r) has two local extrema at 

M 1 a— M 



r 1 = — + -^/M2-2a 1 ^q 2 , r 2 = — - -^Af 2 - 2a X K 2 q 2 . (34) 
We find the number of event horizons depending on the critical values of the parameters a% 



and ro by looking at the behavior of the derivative of the above expression (33) . Thus, n'{r) 



has one root (ri) for a\ < 0, two roots for < a\ < 2^2, one root for a\ = ^^2, {r% = r 2 ) 
and no roots for a\ > 2k2 2 . The number of the horizons will differ according to the sign 
of the function n(r) at these critical points and the values of a\ and r . We will determine 
these intervals and plot the graphs of the function n(r) in the interval for certain values of 
the parameters as follows. 

Case A: For a\ < 0, we can determine the interval of ro related with the numbers of 
horizons from the following inequalities (for a specific case see FIG.l ); 

1. If n(ri) < 0, there are two horizons. 

2. If n(ri) = 0, there is one horizon. 

3. If n(ri) > 0, there is no horizon. 
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FIG. 1: The graph of the function n{r) for the Case:A with a\ = —1, A-1) r$ = 1, 
n[r\) ~ —2.17, A-2) tq ~ 8.83, n(r*i) = 0, A-3) ro = 100, n(n) ~ 2.43 and Reisner-Nordstrom 
case is plotted by the dotted curve (k =M=q=l). 



1. If i) n(ri), n{r2) < or ii) n{ri),n{r2) > there is one horizon. 

2. If i) n(ri) = and nir^) > or ii) n(ri) < and nir^) = there are two 
horizons. 

3. If i) n(ri) < and n(r 2 ) > or ii) n{r\) > and n(r 2 ) < there are three 
horizons. But it is not possible to find an r in the second interval n(ri) > and 
n(r 2 ) < 0. 

Case C: For 1) a>i = 2 ^{ 2 2 and 2) a\ > 2 ^ 2 2 (see FIG. 3); n{r) increases monotonically from 
— oo to oo and there is one horizon. 

The case with no horizon exhibits a naked singularity. The cases with a single horizon 
resemble the extreme Reissner-Nordstrom solution, while the cases with two horizons corre- 
spond to the the generic Reissner-Nordstrom geometry. The case with three horizons seems 
to be new. 

Further investigation of the black hole properties can be done numerically for each of the 
cases above. We plan to deal with this problem in a future work. 
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B-2-i 
B-2-ii 

Reissner-Nordstrom 



(a) for 01 = 0.25 



(b) for ax = 0.25 




Reissner-Notdst 



(c) for a x = 0.25, r = 0.4 



FIG. 2: The graphics of n(r) for the Case-B, B-l-i) ro = 1, n{r\) ~ —0.45, nfa) — —0.19, B-l-ii) 
r = 0.1, n(n) ~ 0.12, n(r 2 ) ~ 0.38, B-2-i) r ~ 0.16, n(n) = 0, n(r 2 ) ~ 0.26, B-2-ii) r ~ 0.47, 
re(ri) ~ —0.26, n(r 2 ) = 0, B-3-i) n(ri) ~ —0.22 and nirz) = 0.04 and Reisner-Nordstrom case is 

plotted by the dotted curve (k =M=q=l). 



IV. CONCLUSION 

We considered a non-minimally y(i?)F 2 -coupled Einstein-Maxwell theory and looked for 
static, spherically symmetric solutions for a specific function Y(R) = 1 ^_ that shows 

1 ^0 

charge screening effects. We obtained a class of asymptotically flat solutions that include new 
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C-1 

C-2 

* * 5 * * * Heissner-Nofdstfem 

FIG. 3: The graph of n(r) C-1) a x = 0.5, r = 1, n(r x ) = n(r 2 ) ~ -0.22, C-2) oi = l,r = 1 and 
Reisner-Nordstrom case is plotted by the dotted curve (« =M=q=l). 



black hole candidate configurations, except for the parameter values a% < and n(r{) > 
when there is a naked essential singularity at the origin. These particular solutions may shed 
light on some problems of gravity such as dark matter and dark energy without introducing 
a cosmological constant or any other exotic fields. This means that, if dark matter is not 
some strange matter, but, for instance the non-minimal couplings produce such effects [23J; 
then the electromagnetic potentials get modified at large (astrophysical) scales and thus 
contribute to the conventional electromagnetic energy density which may then be interpreted 
as the effects of dark matter. We also note that the electric charge q need not have a large 
value to have observable effects. If q is small, a% may be large so that the product aiq 2 
becomes important at large (astrophysical) scales. This can explain the rotation curves of 
galaxies for certain parameter values [23J. 
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